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ABSTRACT 

Recent numerical results on the fractal structure of two-dimensional quan- 
tum gravity coupled to c = — 2 matter are reviewed. Analytic derivation of the 
fractal dimensions based on the Liouville theory and diffusion equation is also 
discussed. Excellent agreements between the numerical and theoretical results 
are obtained. Some problems on the non-universal nature of the fractal structure 
in the continuum limit are pointed out. 



Talk given at 7th Nishinomiya- Yukawa Memorial Symposium, held at 
Nishinomiya-city on 18-20, Nov. 1992, to be published in the Proceedings. 



1. Intoduction 



One of the most important recent developments in gravity theory is that 
we obtained a well defined regularization of quantum gravity in two dimensions. 
This recognition comes from the fact that the continuum formulation and the 
dynamical triangulation are equivalent. The dynamical triangulation can be 
identified as a lattice regularization of quantum gravity and can be analyzed ana- 
lytically by the matrix model and numerically by the computer simulation. ^^'^^ 
It has become clear numerically that the fractal structure of space time is the 
fundamental nature of the quantum gravity. ^^'^^ Analytical investigations by Li- 
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ouville quantum gravity support the numerical results. ' ' ' It is an interesting 
question if the dynamical triangulation works clS cl TG gularization scheme of three- 
and four-dimensional quantum gravity, which may be tested only by numerical 
simulations at this moment. 

In this manuscript we will summarize the recent numerical simulations of two- 
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dimensional quantum gravity coupled to c = —2 matter. ' We also show that 

f8l 

our analytical investigations by Liouville theory have excellent agreements with 
the numerical estimation of a fractal dimension and the mean squared distance 
of gravitational random walks treated by a diffusion equation. 

Our recent analytic investigation on the fractal structure of c = model 
suggests that we must be careful about non-universal behaviors of the fractal 
structures in the continuum limit of dynamical triangulation. 



2. The shape of Typical Quantum Gravitational Space-Time and 
the Fractal Structure 

In the dynamical triangulation of two-dimensional gravity, the metric in- 
tegration in the sense of path integral is replaced by a summation of all the 
different types of surface configuration with a given number of triangles and a 
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given topology. Here we specify the space-time topology as a sphere. Formally 
the continuum partition function is given by 



where Zm[g] is a matter part of the partition function with gravitational back- 
ground and A is the total area. Reguralized counterpart of the above partition 
function by dynamical triangulation is 



where N is the number of equilateral triangles and is the area of the triangle. 
G denotes a triangulation and Gq is the typical triangulation which we select 
from the huge set of triangulations. The last approximate equality in Eq.(2) is 
valid up to the normalization factor and if the selection of the typical surface 
is carried out by a correct procedure which we explain later. Since the path 
integration of the metric is carried out after the selection of the typical surface, 
Go carries the information of the quantum fluctuation of space time effectively. 

A natural question now is: "How does the typical surface look like ?" Since 
the typical surface carries the information of quantum gravitational fluctuation, 

it may look quite different from the classical space time. Some of the possible 
surfaces wc can imagine beforehand may be; 1) smooth surface, 2) spiky surface, 
3) branched surface, .... The next question is: "How do we parametrize the shape 
of the typical surface ?" 

Suppose we obtain the typical surface by the recursive sampling method 
which we explain later, the actual measurements for the typical surface are car- 
ried out as follows. We first fix a marking site on the triangulated surface, 
from which we measure the internal geodesic distance r. We count the num- 
ber of triangles {Vt{r)), links (V/(r)), and sites {Vs{r)) within r geodesic steps. 




(1) 





(2) 



G 
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Vt{r) is essentially the two-dimensional volume in radius r. The number of the 
disconnected boundaries A^6(r) can be counted by using the following relation: 
X = 2 — A^ft(r) and x = Vs{r) — Vi{r) + Vf(r), where x is the Euler number of 
sphere with Nf){r) boundary holes, i^b(r) = Vt{r + 1) — Vt{r) is essentially a 
derivative of Vt{r) with respect to the geodesic distance and is the total length 
of the iVft(r) disconnected boundaries at the distance r from the marking site. 

In the actual calculations we define the geodesic distance on the triangulated 
lattice. To be precise, the region of geodesic distance r from a site is defined as a 
thin connected area composed of all the triangles that are attached to the links 
located at r-steps from the marking site. The number of triangles of this area 
and the number of disconnected boundaries precisely correspond to Li,{r) and 
Ni,{r) respectively. It is also possible to define the geodesic distance on the dual 
lattice. 

We parametrize the above mentioned three quantities as follows: 
{A) Ni,{r) = (number of boundaries at the step r) ~ r", 

(B) Li,(r) = (total length of boundaries at the step r) ~ r^, 

(C) Vt{r) = A{r) = (number of triangles within r steps) ~ r^. 

a, P, and 7 may be called the fractal dimensions of internal space-time geometry if 
they become constant in the large r asymptotic region. If the surface is smooth, 
these quantities behave as follows: A^&(r) ~ 1, Lf,{r) ~ r, Vt{r) ~ r^. Thus 
deviations from these behaviors signal the fractal nature of quantum gravitational 
fluctuation of space-time. In other words we expect that the space-time surface 
is smooth classically while the fractal nature of the space time is an essential 
feature of quantum gravity. 
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3. c = -2 Model and Recursive Sampling Method 

In the numerical simulation of dynamical triangulation there are essentially 
two methods to generate the typical surface of two-dimensional quantum gravity 
: 1) flipping method by Monte Carlo and 2) recursive sampling method which 
is first proposed by Agishtein and Migdal'^^ for pure gravity (c = model). 
The second method can generate much larger number of triangles than the first 
method. The second method, however, necessitates an analytic formula to gen- 
erate a typical surface and is restricted to particular models such as c = and 
c = — 2 models. Here we investigate c = — 2 model which has much simpler 
analytic formula than c = model and is thus easier to simulate larger number 
of triangles: number of triangles = 5 x 10^ for c = —2 while 1.3 x 10^ for c — 0. 

The c = — 2 model or equivalently the two-dimensional gravity coupled to 
c = — 2 matter was introduced and solved analytically in [12]. The partition 
function of two-dimensional gravity coupled to c = — 2 matter is given by 

G <ij> 

= 5]detA(G) (3) 

G 

^ AT + 2'^"^"^^^^/^+^' 

where and t/^^ are fermion fields sitting at a site i of dynamically triangulated 
lattice G. Here we fix the topology of the surface as sphere and the number of 
triangles as N and the area of the triangle is unity. The fermion integration 
leads to a random lattice version of Laplacian A(G). In this particular model 
the determinant of the Laplacian and dual planar 0^ diagrams are related. Tjv+i 
is the number of rooted dual tree diagrams with N + 2 external legs while Rn /2+i 
is a number of rainbow diagrams with N/2 + 1 lines. A rainbow diagram can be 
constructed by connecting the A'" -|- 2 external legs of a tree diagram in such a way 
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that none of the N/2 + 1 hnes crosses over. The third equahty in eq. (3) is a very- 
useful one and can be estabhshed by an apphcation of the well-known Kirchhoff 
theorem. The introduction of the fermion field effectively means an embedding 
of a surface into c = — 2 dimension. 

Tn satisfies the Schwinger-Dyson equation: Tn = Ylk=i '^kTn-ki which can 
be solved by introducing a corresponding generating function of T„ and leads to 
the solution: T„ = (2n — 2)!/n! (n — 1)!. Rn satisfies the same type of relation 
as Tn, which leads to a relation: Rn-i = T^. 

The recursive sampling algorithm to generate a typical surface goes as follows. 
First of all one of the great advantages of c = —2 model is the factorization 
feature of Tjv+i and Rn/2+i i^i Eq.(3). In order to generate a typical surface 
we independently generate a tree diagram and a rainbow diagram with a correct 
probability and then connect them. The branching probability to divide a tree 
diagram with n + 1 legs into two different tree diagrams with k-\-l external legs 
and n — k + 1 external legs is given by W{n, k) = T^Tn-k/Tn- A typical tree 
diagram with N vertices can be generated by applying N times of the probability 
formula W{n, k). A typical rainbow diagram can be generated similarly. 

4. Numerical Results on Fractal Structures of c = -2 Model 

We show the numerical results of the fractal dimensions a, [3 and 7 with 
the fittings of (A) Ni,{r) ~ r", (B) Li,{r) ~ r^, and (C) Vt{r) = A{r) ~ r'^ 
in Figs.l, 2, and 3, respectively. We have measured those quantities for the 
following number of triangles: (1) 8 x 10^, (2) 4 x 10^, (3) 2 x 10^ (4) 10^, and 
(5) 5 X 10^. The values a, /? and 7 are approaching to constant values 2.55 ± 0.1, 
2.53 ± 0.1 and ~ 3.5, respectively, and thus show clear fractal behaviors. We 
obtain an approximate relation 7 — /3 ~ 1 numerically, which is expected from 
the relation dVt{r)/dr ~ Li){r). What is interesting and unexpected is that the 
two-dimensional quantum gravitational space time is very much branching. For 
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example the space-time surface with c = — 2 matter sphts into ~6000 branches 
after ~60 geodesic steps in the case of 5 x 10^ triangles. See Fig. 1(a). 

It happens that the fractal dimension of the number of boundaries (a) is 
approximately same as that of the total length of the boundaries To es- 

timate the average perimeter length of a branch we have measured the ratio 
(D): Li,{r) / Nf,{r) in Fig.4. It shows clear constant behavior as a function of the 
geodesic distance. The numerical value is approximately ~15 in the case of 10^ 
triangles. It looks as if all the branches have similar perimeter length. In order 
to see if this is the case we have measured (E): the number of boundaries Pr{l) 
with a given boundary length / as a function of the boundary length / for a given 
geodesic steps r measured from a marking point. The results are shown in Fig. 5. 
As we can see in the figures, different sizes of boundary length are uniformly dis- 
tributed with roughly a power law behavior of ~ In other words there is no 
particular preference of the perimeter length even though the average perimeter 
length is ~ 15. 

We can now imagine how the two-dimensional quantum gravitational space 
time looks like. It has many branches whose perimeter length varies small size to 
large size. Number of the branches with large perimeter size is small but there 
exists such branches with a certain probability. If we look into this quantum 
gravitational universe microscopically, we feel like being in a dense jungle and 
the universe looks like complete chaos. The universe, however, have beautiful 
fractal structures and thus has an order as a whole. We have given here a 
numerical results of c = —2 model. In c = model Agishtein and Migdal found 
a similar branching behavior of the space-time surface, but didn't observe the 
fractal behavior. This suggests that the branching behavior of the quantum 
gravitational space time is quite universal phenomena irrespective of the choice 
of matter. 

The geodesic distance so far has been measured on the triangulated lattice. 
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Geodesic distance on the dual lattice or equivalently on the dual planar (f)^ di- 
agram is defined as a number of minimum steps located between two vertices. 
It is one of interesting questions whether the numerical results of the fractal 
structure change depending on the different definitions of the geodesic distance. 
We show the numerical result of the fractal dimension of boundary length (F): 
/3(r) = (i log L;,(r)/(i log r measured with the geodesic distance on the dual lattice 
in Fig. 6, which should be compared with the corresponding result in Fig. 2. As 
we can see, the behavior to approach the stationary value is very slow in the case 
of geodesic distance on the dual lattice. This may be understood by recognizing 
the fact that the geodesic distance on the dual lattice generally takes much more 
steps than those on the triangulated lattice to measure the same distance. It 
is, however, natural to expect that P defined on the dual lattice will eventually 
approach to the same value as that defined on the triangulated lattice. 

5. Diffusion Equation on the Dynamically Triangulated Lattice and 
the Continuum Limit 

In this section we try to investigate the fractal structure of quantum gravity 
analytically by diffusion equation. We first define adjacency matrix Kij. For a 
given typical surface Gq we number the sites of the triangulated lattice. Then the 
(i, j) component of the adjacency matrix Kij is one if i-th site and j-th site are 
connected by a link as next neighboring sites, and zero if they are not connected 
by a single link. It is interesting to note that (n, uq) component of {K^)n,no 
counts the number of possible random walks reaching from a marking site uq to 
a site n after T steps. The Laplacian defined on the dynamically triangulated 
lattice is given by 

Al ^ 1- S, Sij = -Kij, (4) 

Qi 

where qi is called coordination number and denotes a number of links connected 
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to the site i. Sij is thus a probabihty of one step random walk from the site j 
to the neighboring site i. The diffusion equation on a given surface Go with N 
triangles is now given by 

ar*S?"^(T;n,no) = Az.(Go) vE'SJ'^^lT; n, no), (5) 

where Ot is a difference operator in T direction and " (T; n, no) is a wave 
function of the diffusion equation and denotes the probability of finding the ran- 
dom walker at the site n after T steps from the starting site no. A solution of the 
diflFusion equation can be easily obtained as j^"-* (T; n, no) = e^^^^'^°\5n,no)i 
where (5n,no) is A'"-component vector with unit no entry. 

We now consider the continuum limit of this diffusion equation. First of all 
we recover the lattice constant a. In taking continuum limit, the total physical 
area A = afNi is fixed and ai —>■ {Ni oo) is taken, where Ni is the number 
of triangles and of is the area of a triangle. In each step of the limiting process 
we select a typical surface Gi for the given number of triangles N^, on which 
the lattice Laplacian A^(Gi) of Eq.(4) is defined. Now the lattice version of the 
diffusion equation (5) can be rewritten as 

^{^f\T + alx,xo) - ^f\T-x,xo)} = ^AL{Gi)^f'\T-x,xo)}, (6) 

where the location of the site x is measured with respect to the lattice constant 
ai. Thus we identify the dimension of T as that of area: dim[T] = dim[A]. In the 
continuum limit the solution of the diffusion equation (6) is expected to approach 
the continuum wave function: '^^^'\T;x,xo) '^f^\T;x,XQ). Numerically 
we approximate the limiting surface as the typical surface (Go) of the maximum 
size: Goo — Gq As we have already noted in Eq.(2), the metric integration is 
effectively carried out for the equation (6) since we have chosen a typical surface. 
This means that the quantum effect is included for the wave function of Eq.(6). 
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On the other hand the solution of the continuum counterpart of the diflFusion 
equation: dr^{T; x, xq) = A{g)^{T;x,xo) is still background metric dependent 
in general. Furthermore the dimensions of T and r may not necessarily be equal. 

5. Fractal Dimensions by Liouville Theory 

An analytical treatment of fractal dimensions by Liouville theory has been 
first given by Kawai and Ninomiya. In their treatment fermion is introduced 
as a test particle to derive the fractal dimension. It has been recognized that 
there are several ways of defining fractal dimensions, which may explain the 
discrepancy between the theoretical and numerical results. Here we derive 
fractal dimensions by investigating the gravitational random walks with the help 
of diffusion equation and Liouville theory. In this section we briefly sketch the 
derivation of the fractal dimension while the details of the derivation will be given 
elsewhere. 

Let us now define the comeback probability of random walk on the triangu- 
lated lattice and relate it with the continuum expression of Liouville theory as 
follows: 



where < O >a is an expectation value of O with the partition function given by 
Eq.(l). We should remind of the fact that the metric integration is effectively 
carried out since we have chosen the typical surface Gq for the wave function of 
the comeback probability. The initial wave function can be formally written as 
\l/(0; ,x', x) = liuiy^x S{y — x)l/ and need to be reguralizcd. The first similarity 
relation of Eq.(7) can be understood from the relation: f dy^Jg'^{T;y,x) = 1, 
which is obvious from the provability interpretation of the wave function. The 
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G(T) 



*^"°^(T;a;o,a;o) 



< J dx^ >A 



(7) 




dimensional arguments coming from the equation (6) supports the last similarity 
relation. Here comes a prediction of comeback provability (G): G{T)T ~ const.. 

We next consider how to accommodate the Weyl invariance into the diffusion 
equation of random walk by using DDK formulation of Liouville theory. Let 
us consider the following quantity by Liouville theory: 

< dx^-^{T;x,x) >A= < dxy/g-^{0;x,x) >A + 

\ (8) 

where the solution of the diffusion equation is expanded by r. 

Taking a conformal gauge g^v{x) = Qf^iyC^^^^ and introducing DDK arguments, 
we can rewrite the first and second terms of the Eq.(8) as 

< J dxy/g'^{0;x,x) >a = < J dx^/^[^S{x - xq)]^^^^ >a, 

< J dx^A^{0-x,x) >A ^< J dx^g [A^ e""^"^ - xq)]^^^^ >a, 

(9) 

where the term e^-i*^ is introduced to keep the Weyl invariance of the second 
term. The expectation value < 0{g) >a is now rewritten by using the well- 
known expression 



< 



0{g) >A= I Vg(j>Z^MZM[g]5{ I dxyfge'^^^ - A) 0{g,cf>) 

D-25 
exp( Si,[4>,g]), 

where Zpp is the Fadeev Popov contribution and Si^[4>,g] is the Liouville action. 
an appeared in Eqs.(9) and (10) is given by 

2n 

Q^n = , (II) 

1 + V(25-c-24n)/(25-c) ^ ^ 

Invariance of the expectation value under the translation of the conformal field 
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(j) ^ (j) — InX/ai leads to the change: 

5(^J dx^ e"!"^ - Aj ^ X5(^J dxy/g e"i<^ - A^^ , which is then interpreted as the 
scale change of the physical area A — > A^. To require the invar iance of each term 
of Eqs.(8) and (9) under this scale transformation r should scale as r — > A "i r. 
In other words, 

"—1 

dimr = dim^ (12) 

We now point out that the expectation value of the mean squared geodesic 
distance is evaluated by the standard continuum treatment 

j dx^{r{x,XQ)f-^{Q;x,x) = -4t + ©(r^), (13) 

which is now related with the quantum version of the mean squared geodesic 
distance in the small r region 

<r2> = Y,{r{x,XQ)f-^^^''\T-x,XQ) 

X 

^ < f dx^ J dxoy/g{r{x,xo)}'^'^{T;x,x) >a (14) 
~ < / dxy/g >A 



where the geodesic distance on the dynamically triangulated lattice is same as 

that of the numerical simulation. The last three similarity relations are due 

to the dimensional arguments. We thus obtain an analytic prediction of the T 

2 

dependence of the mean squared geodesic distance (H): < >~ T «i = Tt(=) . 

From the dimensional arguments of Eq.(14), we obtain the following relation: 

dimlA] = (izm[r~^"i/"-i] = dimlr"'^^^]. In the numerical simulation (C): l^l'^) = 
A{r) ~ r^, the fractal dimension 7 for c = —2 is estimated from the measurement 
of the two-dimensional volume of space time. It is then natural to expect that 
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the dimension of A{r) and that of the total area A is same. We then obtain the 
analytic evaluation of the fractal dimension 



f\ o "1 n.. V25 - c + V49 - c 

a-i V25 - c + V 1 - c 

where some of the typical values of the fractal dimension are; 7(1) = 
2(1 + v^), 7(0) = 4, 7(-2) = i(3 + Vl7) = 3.56 7(-oo) = 2. We show 
the c dependence of the fractal dimension 7(c) in Fig. 7. As we can see from 
Fig. 7, the fractal dimension varies smoothly from the classical value 7(— cxd) = 2 
to 7(1) = 2(1 + \/2) and gets imaginary for c > 1. 

We now summarize the analytic predictions: 

(G) G{T)T ~ const., 

(H) < > = T A ~ 2^-0.56^ 

(/) 7(-2) = i(3 + 717) = 3.56 

First of all the analytic prediction of 7(— 2) = 3.56 ■ ■ ■ should be compared with 
the numerical value of 7 ~ 3.5. The theoretical value (0.561) has an excellent 
agreement with the experimental value. We show the numerical evaluations of 
(G): G{T)T and the mean squared distance of random walks with gravitational 
quantum fluctuations (H): < > in Fig.8 and Fig.9, respectively. As we can 
see, G{T)T in Fig.8 show the clear constancy and thus excellent agreement with 
the prediction (G). In Fig.9 the numerical value of the power of T is approaching 
to the theoretical value and slightly away in the large T region where finite size 
effects may be important. If wc consider the accuracy of the vertical measure 
we may conclude that the agreement with the theoretical prediction is excellent 
again. 
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6. Non-universal Nature of the Fractal Structure in the Continuum 
Limit 

In our recent investigation '^^^ we have obtained analytic formulation to eval- 
uate the fractal structure of two dimensional quantum gravity without matter 
(c = 0). In this analysis we have used the results of matrix model and introduced 
a combinatorial consideration and then succeeded to derive a transfer matrix. 
The notion of the geodesic distance explained in this manuscript played an im- 
portant role in the formulation. 

We briefly summarize the main conclusion of this investigation. ^^^^ Let us 
first define a continuum function p(L; D) which is a function of the boundary 
length L of a loop located at the continuum geodesic distance D measured from 
a marking point. p{L;D)dL counts the number of boundaries whose boundary 
lengths lie between L and L + dL. It is evaluated by taking the continuum limit 
from the transfer matrix and disk amplitude of dynamical triangulation. The 
functional form of p{L; D) for c = model is given by 



(b) arc the lattice counterparts of the multiplicity function of c = —2 model; 
(a) p{L,r = 18) and (b) p{L,r = 40), respectively, where continuum geodesic 
distance D is replaced by the lattice geodesic step r. In Fig. 5 only the small x 
region is shown. 

In order to examine the scaling property of the fractal structure, it is conve- 



P{L; D)D^ = + 1-3/2 + ^^i/2)e-, 



(16) 
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nient to introduce the following quantities: 

/■oo 

<L">= / dLU'p{L-D). (17) 

We can then derive the fractal scaling behaviors of c = model on the following 
quantities 

< L° > = Ni,{D) ~ const x D^a'^/"^ + const x Da~'^l'^ + const x 

< > = Lb(i:>) ~ const X D^a''^/'^ + const x (18) 

< L'' > - const x L>2« > 2), 

where a is the lattice constant and Ni,{D) and L}){D) are the same quantities as 
those defined in (A): A^ft(r) ~ and (B): Li,{r) ~ r''^ in the numerical simulations 
except that the argument is now the continuum geodesic distance D. 

As we can see from Eq.(18), A^b(-D) and Lf){D) include the inverse power of 
the lattice constant dependent part as a dominant contribution. For this c = 
model the fractal dimension 7 = /5 + 1 = 4 obtained from the D dependence of 
Li,{D) coincides with that of 7(0) = 4 obtained from the formula (15). What is 
unexpected is that the coefficient of this power dependent term is lattice constant 
dependent and thus has non-universal nature. It should also be noted that the 
fractal dimension obtained from < >~ D^, which has the same dimension as 
A{D) ~ D'^, happens to reproduce the same fractal dimension 7 = 4 discussed 
above. This fractal dimension derived from the < > does not have lattice 
constant dependence and thus should have a universal nature. It is a crucial 
question if 7(0) = 4 obtained from the Liouville theory is equivalent either with 
the one derived from < >= Li,{D) or with another one derived from < >. 

The analytic result of c = model given in Eq.(18) suggests that A^b(-D) 
and Li,{D) include non-universal part as a dominant contribution which become 
irrelevant in the continuum limit but show a fractal nature in the microscopic 
level. 
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Agishtein and Migdal have carried out the numerical simulation for c = 
model and measured the quantity ^{r) = d log A(r) j d log r with the maximum 
number of triangles ~ 1.3 x 10^. The observed fractal dimension at the maximum 
number of traingles is roughly 3 and still far below the analytically expected 
value 4. Wc believe that 7(r) will approach to the analytic value for much larger 
number of triangles. 

7. Conclusion and Discussions 

In this manuscript we have reported the following numerical results and an- 
alytic predictions and corresponding figures: 

(A) Ni,{r) = (number of boundaries at the step r) ~ r" Fig.l, 

(B) Lf,(r) = (total length of boundaries at the step r) ~ Fig,2, 

(C) Vt{r) = A{r) = (number of triangles within r steps) ~ Fig. 3, 

(D) Lbir)/Nbir) Fig.4, 

(E) The number of boundaries Pr{l) with a given boundary length / as a func- 
tion of the boundary length I for a given geodesic step r measured from a 
marking point Fig. 5, 

(F) Lh{r) ~ r^, where r is defined on the dual lattice Fig, 6, 

(G) G{T)T ~ const. -Fig.8, 

(H) < > = T A ~ -Fig.9, 

(/) 7(-2) = i(3 + Vrf) = 3.56 -Fig.3. 

The numerical results (A), (B) and (C) show clear fractal structure of two- 
dimensional quantum gravitational space-time. In particular the space time is 
violently branching. The results (D) and (E) show that the perimeter length of 
the branches varies from small to large sizes but the average size is independent 
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of the geodesic distance measured from a marking point. This suggests an ex- 
istence of some rule even for c = — 2 model based on the fractal nature of the 
branching behavior . In fact for pure gravity (c = 0) we have analytically de- 
rived the c = counterpart of the multiplicity function (E) of Fig. 5. The result 
(F) in Fig. 6 shows that the geodesic distance defined of the dual lattice enforces 
the very slow approach of the fractal dimension to the asymptotic value. The 
analytic results (G) and (H) have excellent agreement with the numerical results 
of corresponding figures Fig. 8 and Fig. 9, respectively. Thus the two-dimensional 
quantum gravity with dynamical triangulation can be treated by the diffusion 
equation of random walk with quantum gravitational fluctuations. The lattice 
version of the wave function and the continuum counterpart of diffusion equa- 
tion are related by Eq.(7), which makes it possible to accommodate the Liouville 
theory. The fractal dimension (I): 7(— 2) obtained from the formula (15) excel- 
lently agrees with the numerical results of c = —2 model obtained from Fig. 2 
and Fig. 3 while the 7(0) obtained from (15) coincides with the analytic results 
derived from (18). Therefore we tend to believe that the formula (15) derived 
from the Liouville theory provides correct fractal dimensions of two-dimensional 
quantum gravity with matter fields. 

As we have mentioned in the last section, Ni,{D) and Li,{D) of c = model 
include non-universal lattice constant dependent part as a dominant power con- 
tribution of geodesic distance which becomes irrelevant in the continuum limit. 
This may also be the case for c = — 2 model. If Nf,{D) and Lf,{D) include non- 
universal part even for c = —2 model, we are puzzled why the fractal dimension 
7 ~ /3 -h 1 ~ 3.5 obtained from the numerical results of Li,{r) and A{r) ex- 
cellently agrees with the analytic result 7(— 2) = 3.56 • • • obtained from the 
formula (15) which correctly reproduces the c = analytic result. This kind of 
question may be partially answered by numerical measurements for the models 
other than c — 0. Recently we have measured p{L; D) and < > for c = —2 
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model numerically and observed the similar scaling behavior as c = model. 
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The problems are not yet completely settled at this moment since we have not 
yet obtained the analytic results of < L"^ > for the model other than c = 0. 
There are still several other issues to be cleared up: in particular the relation 
between the derivations of fractal dimensions by Liouville theory and the exact 
treatment. 
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FIGURE CAPTIONS 



Fig. la (A): A^fe(r) — r dependence, where N}j{r) is the number of boundaries at 
the step r for various number of triangles of sphere topology: (1) 8 x 10^, 
(2) 4 X 10^, (3) 2 X 10^ (4) 10^ and (5) 5 x 10^ 

Fig. lb (A): a{r) = d\ogNi,{r) / dlogr for various number of triangles, where a{r) 
is the fractal dimension parametrized in (A). 

Fig. 2 (B): (3{r) = dlogLf){r)/dlogr for various sizes of triangles, where the 
fractal dimension /3(r) is parametrized in (B). 

Fig. 3 (C):7(r) = d\ogVt{r)/dlogr for various sizes of triangles, where 7(r) is 
parametrized in (C). 

Fig. 4 (D): Lf,(r)/Ni,(r) for the case of 10^ triangles. 

Fig. 5 (E): The number of boundaries P{1) with a given boundary length / as a 
function of the boundary length I for the geodesic distance (a) r = 18 and 
(b) r = e = 40, where the number of triangles is 10^. 

Fig. 6 (F): /3{r) = dlog Li,{r) / dlogr for various sizes of triangles, where the 
geodesic distance here is defined on the dual lattice. 

Fig. 7 The fractal dimension 7(c) given by Eq.(15) as a function of the matter 
central charge c. 

Fig. 8 (G): G{T)T as a function of geodesic step T, where G{T) is the comeback 
probability of random walk given by Eq.(7) for the case of 10^ triangles. 

Fig. 9 (H): 5 = dlog < > /dlogT as a function of geodesic step T, where 
< > is the mean squared geodesic distance. The theoretical value 0.561 
is shown as a solid line in the figure. 
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